Fermi Surface and Order Parameter Driven Vortex Lattice Structure Transitions in 

Twin-Free YBaaCusOy 
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We report on small-angle neutron scattering studies of the intrinsic vortex lattice (VL) structure 
in detwinned YBa2Cu307 at 2 K, and in fields up to 10.8 T. Due to the suppressed pinning to twin- 
domain boundaries, a new distorted hexagonal VL structure phase is stabilized at intermediate 
fields. It is separated from a low-field hexagonal phase of different orientation and distortion by a 
first-order transition at 2.0(2) T that is probably driven by Fermi surface effects. We argue that 
another first-order transition at 6.7(2) T, into a rhombic structure with a distortion of opposite sign, 
marks a crossover from a regime where Fermi surface anisotropy is dominant, to one where the VL 
structure and distortion is controlled by the order parameter anisotropy. 

PACS numbers: 74.25.Qt, 74.72.Bk, 61.05.fg 



YBa2Cu307_5 (YBCO) continues to be of great fun- 
damental interest for research into High-Tc materials. 
In particular, the study of the m^netic vortex lattice 
(VL) commands attention [H, S S HI , as the VL struc- 
ture and distortion directly reflect the underlying mi- 
croscopic state and inter- vortex interaction. In YBCO, 
superconductivity on the CuO chains, which lie along 
the crystal b-axis, is believed to be induced by proxim- 
ity effect from the Cu02 planes [1, 0]. This results in 
an enhanced supercurrent response in this direction rel- 
ative to the a-axis, which is reflected in the distortion 
of hexagonal VL structures seen at low fields [2, 3, 4]. 
At low magnetic fields parallel to the c-axis {B \\ c), 
anisotropic London theory relates the value of this dis- 
tortion directly to the penetration depth anisotropy jab 
(= Xa/Xb = \/^a/^b) B Isl- However, other theories 
also predict that the VL structure and distortion depend 
on the effects of a nodal gap structure [10, 111. [12I. [13I . 
and on non-local electrodynamics combined with a Fermi 
surface anisotropy [15]. The common prediction is for the 
formation of square VL structures at high fields. How- 
ever, for YBCO specifically, comparison to theory is com- 
plicated by the orthorhombic crystal structure and in- 
plane electronic anisotropy. These will introduce an 5- 
wave admixture in phase with the predominantly dj.2_y2 
order parameter, th oug h the exact details remain under 
debate [l^, 17, [3, fioj . Furthermore, previous studies 
of the VL in YBCO have been hampered by twin plane 
pinning along {110} [H, HI, 0, 0, HI , which may obscure the 
intrinsic VL structure. For the measurements reported 
here, the sample was detwinned and overdoped {5 = 0) 



to minimize vortex pinning. As a result, for the first time 
we have observed VL structures, with || c, that avoid 
the effects of twin plane pinning. 

To image the VL structure directly, we have used the 
bulk probe of small-angle neutron scattering (SANS) in 
a series of experiments at the ILL and PSI. Neutrons of 
wavelengths between 6 - 10 A were selected with a 10% 
FWHM spread and collimated over 6 - 14 m before the 
sample. Diffracted neutrons were collected by a posi- 
tion adjustable 2D multi-detector. The sample, of total 
mass ~ 30 mg, was a mosaic of six single crystals, co- 
aligned to ^ 1.5° about the c-axis. The crystals were 
grown in a BaZrOs crucible [20] , detwinned under uniax- 
ial stress at ^ 500°C for 24 hours [21], and oxygenated 
to O7 under high pressure [22[. Using neutron diffraction 
to observe (100) and (010) reflections from the sample, 
no minority twin-domain signal was detected at the 1 % 
level. The mosaic was mounted with the a-axis verti- 
cal and the c-axis parallel to the neutron beam. With 
5||c, the VL was prepared using two different techniques. 
The first was a standard field-cool through Tc in a con- 
stant amplitude field ("field-cooling", FC). The second 
was a field-cool in a field oscillating around the target 
value ("oscillation field-cooling", OFC) with an ampli- 
tude between 0.1 % and 0.2 % of the target field, and a 
frequency oi ^ 2 min~^. For both FC and OFC, once 
at the measuring temperature of 2 K, the field was held 
stationary for diffraction measurements. These were per- 
formed by rotating the cryomagnet and sample together 
to angles that brought various diffraction spots onto the 
Bragg condition at the detector. In all cases, background 
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FIG. 1: (Color online). The VL diffraction patterns taken at 2 K in applied fields, B of (a) 1.5 T after FC, (b) 1.5 T after 
OFC, (c) 4 T after FC, (d) 6.5 T after FC and (e) 10.8 T after FC. Overlaid patterns indicate the different VL structures that 
make up the overall diffraction patterns, and the angles between certain Bragg spots. In (b) the angle is bisected by the b* 
direction, in (c) p is bisected by a* whilst in (d) v is bisected by b*. In all cases, white arrows indicate the {110} directions. 
The diffraction patterns were constructed by summing the detector measurements taken for a series of sample angles about 
the horizontal and vertical axes. Statistical noise around the main beam at the center of the image has been masked, and the 
pixelated data has been smoothed with a Gaussian envelope smaller than the instrument resolution. The real space VL can be 
visualized by rotating the reciprocal space image by 90° about the field axis and adding an additional spot at the center. 



measurements were taken at T>Tc, and subtracted from 
2 K foregrounds, to leave just the VL signal. 

Figs. 1(a) to (e) show VL diffraction patterns col- 
lected at 2 K and at various fields. The results at 1.5 T 
(Fig. 1(a) and (b)), reveal the marked difference in VL 
structure we observe depending on whether the VL at low 
fields is prepared by FC or OFC. In Fig. 1(a), close in- 
spection of the VL structure reveals it to be composed of 
two distorted hexagonal domains, each with two of their 
six spots corresponding to {110} planes. These struc- 
tures are indicated in Fig. 1(a). Similar VL structures 
arise in one of the crystal domains of a twinned crys- 
tal 0, with the VL orientation controlled by twin planes. 
In the present case, there was no detectable signal from 
VLs with a distortion corresponding to the other crys- 
tal domain. Nevertheless, we cannot rule out pinning 
by twin planes associated with a very small fraction of 
minority crystal domains in our sample. However, the 
pinning potential may instead be due to minor crystal 
disorder along {110} planes remaining after the removal 
of twins. After OFC at 1.5 T, Fig. 1(b) shows that the 
VL structure is instead composed of a single distorted 
hexagonal domain. The structure change is also accom- 
panied by a > 50% mean reduction in the rocking curve 
width of the Bragg spots, and > 20% mean reduction 
in their azimuthal width. These correspond to substan- 
tial improvements in the longitudinal correlation of the 
vortices in the field direction, and in the orientational or- 
der about the field axis respectively. This improvement 
in VL perfection suggests that the vortex motion caused 
by OFC allows the VL to explore the true minimum 
of its free energy F, at temperatures below the FC irre- 
versibility line. Since anisotropic London theory predicts 
all orientations to be degenerate Q, we expect to be 
shallow at low fields, and OFC required to reveal it. 

At fields from 2.5 to 6 T, with stronger vortex- vortex 
interaction, the effect of OFC is observed to be less im- 



portant; whether the VL is prepared by FC or OFC, the 
VL structure is composed of a single distorted hexag- 
onal domain like that in Fig. 1(c). Note that this 
intermediate-field structure is rotated by 90° relative to 
that in Fig. 1(b). However, in the field region between 
1.5 and 2.5 T, the structures illustrated in Figs. 1(b) and 
(c) will have very similar - and even more shallow - F^, 
so the free energy F must vary little with VL orienta- 
tion. Under these circumstances, pinning may again be 
expected to play a role. Indeed, on FC at 2 T, we observe 
that the diffracted intensity forms an almost continuous 
ellipse, indicating VLs pinned in a range of orientations 
of similar F. Even with OFC at 2 T, due to the absence 
of a dominant intrinsic F^, pinning remains important, 
as we observe the re-emergence of the {110}-pinned VL 
structure of Fig. 1(a). However, now this structure co- 
exists up to 2.25 T with that illustrated in Fig. 1(c). In no 
case do we observe a smooth transition between the two 
VL orientations represented by Figs. 1(b) and (c), which 
indicates that the intrinsic transition between them is 
likely to be first- order. Further evidence for this is given 
by the field-dependence of the hexagonal structure dis- 
tortion. For each VL structure, this is characterized by 
the axial ratio, r^, of the ellipse that overlays the Bragg 
spots in each domain. The value of 77 may be calculated 
from the measured angles between VL reciprocal lattice 
vectors (all angles would be 60° for the isotropic case of 
T] = 1). Angles (approximately) bisected by symmetry 
axes (i.e. uj in Fig. 1(a), (j) in Fig. 1(b) and p in Fig. 1(c)) 
provide the most sensitive measures of r]. The resulting 
field-dependence of r] is presented in Fig. 2. We note that 
for the OFC data. Fig. 2 reveals a clear discontinuity in 
?7 at ~ 2 T, which is further evidence for the first-order 
nature of the re-orientation transition. 

Further increase of the field yields another first-order 
VL structure transition. Fig. 1(d) shows that by 6.5 T 
the single distorted hexagonal structure like that in 
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FIG. 2: (Color online). The field- dependence of the axial 
ratio, rj of the ellipse that overlays the Bragg reflections in 
each VL domain. Triangles, diamonds and circles represent 
structure types illustrated in Figs. 1(a), (b) and (c) respec- 
tively. Open and filled symbols represent VLs prepared by FC 
and OFC respectively. Error bars (including those in Fig. 3) 
indicate the standard error. For a VL structure like that in 
Fig. 1(a), we present the mean value of 77 for the two domains. 
The blue line is a guide to the eye for the OFC data for the 
VL structure like that in Fig. 1(c). The shaded portion indi- 
cates the field range in which we propose the first-order VL 
structure transition occurs. 



Fig. 1(c) co-exists with a newly emerged domain, that we 
refer to as 'rhombic' Note that for the distorted hexag- 
onal structures, there are six spots of comparable inten- 
sity, and only four in the rhombic phase. Remarkably, 
Fig. 1(d) also reveals that the long axis of the distortion 
of the rhombic diffraction pattern lies along b*, oppo- 
site to that of the hexagonal structure which lies along 
a*. The clear co-existence of the two structures at 6.5 T, 
and the rapid crossover of diffracted intensity from one 
domain to the other, indicate the transition to be first- 
order. At 7.5 T, the VL in the entire sample is composed 
solely of the rhombic structure, whose aspect ratio falls 
monotonically towards unity with increasing field. By 
10.8 T (Fig. 1(e)), the vortex nearest neighbor directions 
are close to (110), showing a nearly square structure. 

An inspection of Figs. 1(b) - (e) shows that none of 
the intrinsic VL structures have VL planes exactly par- 
allel to {110}. In all previously investigated samples 
(1, 0, 0, m, one plane in each VL domain was pinned 
to a {110} twin plane, manifest or residual, for any field 
B \\ c. As a consequence, those VL structures evolve con- 
tinuously at high field into the nearly square structure, 
maintaining the effects of twin plane pinning [3]. In our 
study, the weakness of {110} pinning allows the intrin- 
sic intermediate-field phase (Fig. 1(c)) to be observed for 
the first time, and to undergo the expected first-order 
transition p^], to a nearly square phase at high field. 

Fig. 3 presents a more detailed view of the VL struc- 
ture and distortion. For each structure type. Fig. 3 shows 
the field-dependence of the apex angle between the basis 
vectors of the primitive cell. The analysis confirms that. 
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FIG. 3: The field-dependence of the VL unit cell apex an- 
gle for the various structure types. For clarity, we only show 
measurements for VL structures aligned with the crystal axes, 
with inset diagrams showing the real space orientation of the 
various VL unit cells. The indicated angles in these diagrams 
correspond exactly to those shown in the diffraction patterns 
of Fig. 1. Diamonds, circles and squares indicate VL struc- 
tures like those in Fig. 1(b), Fig. 1(c) and Fig. 1(e) respec- 
tively. Open and filled symbols represent VLs prepared by 
FC and OFC respectively. Error bars not visible are compa- 
rable with the symbol size. The dotted line indicates an apex 
angle of 90° for a perfectly square structure. The dash-dot 
line indicates an apex angle of 60° for an isotropic hexagonal 
structure. Shaded portions indicate field regions over which 
we observe first-order VL structure transitions. 



in spite of the effect of OFC which appears to promote 
a more square structure at high fields, even at maximum 
field we do not observe a perfectly square VL. It is pos- 
sible that such a structure occurs beyond our available 
field range. However, since the crystal is orthorhombic, 
there is no symmetry reason for the stabilization of a 
square VL structure. Fig. 3 also shows the two shaded 
field regions over which we observe first-order structure 
transitions. The size of these regions probably indicate 
slight sample inhomogeneity, with different structures be- 
ing frozen-in in different parts of the sample. From these 
however, we estimate the low- field transition occurs at 
2.0(2) T, and the high-field transition at 6.7(2) T. Prelim- 
inary measurements indicate that these transition fields 
are not strongly T-dependent. Fig. 3 also shows that 
the rhombic structure phase has more noticeable differ- 
ences in the precise shape of VLs prepared by FC and 
OFC, than in the intermediate field hexagonal structure 
phase. This apparent change in i^-dependence indicates 
that the transition at 6.7 T coincides with a change in 
the mechanism controlling the VL structure alignment. 

As noted earlier, anisotropic London theory for high- 
K, superconductors ^ cannot explain the first-order re- 
orientation transition that we observe at 2 T. More 
promising alternatives lie in extensions to the anisotropic 
London model containing non-local corrections that cou- 
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pie the VL to properties of the Fermi surface. Examples 
are an 5-wave theory [15|, or a similar theory allowing 
for a d-wave gap [23|. Each predicts first-order 45° re- 
orientation transitions when the field is parallel to a four- 
fold symmetry direction. In the two-fold symmetry of 
YBCO, we observe instead a first-order 90° re-orientation 
transition. To overcome the detailed differences between 
prediction and experiment, it is likely that higher-order 
terms are required. 

At higher fields however, the first-order transition at 
6.7 T is more likely to be driven by the (i-wave order pa- 
rameter. Notably, the d-wave theory of predicts the 
transition between hexagonal and square VL structures 
to be first-order, as their free energies cross. Currently 
available theories 13, [HI, [H, [H, [13 treat tetragonal sys- 
tems with a pure (i-wave order-parameter, but we expect 
that YBCO is close enough to this limit that valid com- 
parisons can be made. At high B and low T, such theories 
predict that square structures are stabilized with the vor- 
tex nearest neighbor directions lying along the nodes of 
the gap. For a pure (1^2 _y2 order parameter, this is along 
(110), in agreement with the directions approached at 
high B by the vortex nearest neighbors in the rhombic 
phase. 

Finally, we discuss the detailed shape of the VL struc- 
tures. Within anisotropic London theory, the hexagonal 
VL structure distortion parameter, rj = jab^ the penetra- 
tion depth anisotropy [9]. On referring to Fig. 2, we note 
that our values from the OFC data at low fields agree 
with previous work [3, [1, [3- However, if the structure 
transition at 2 T is driven by non-local effects, this will 
render 77 7^ 7^5, and restrict the local regime only to the 
lowest fields. Nevertheless, above 2 T, the monotonic fall 
of T] with increasing field does suggest the suppression of 
an underlying anisotropy. This could be related to a pos- 
sible field-driven suppression of the proximity-effect in- 
duced superconductivity on the chains [7]. Importantly, 
at 7 T there remains a significant a-b anisotropy that 
is likely to persist to higher fields. Hence, on moving 
through the transition at 6.7 T, the change in the sign 
of the distortion between the hexagonal and rhombic 
structures is unlikely to be described by a Fermi sur- 
face anisotropy. Instead we consider the structure of the 
order parameter. Recent very low field phase-sensitive 
measurements [l8| show that the dx2-y2-\ike nodes do 
not lie exactly along the (110) directions, but at ?^±50° 
about the b-axis. Such an imbalance in the lobes may be 
described by a d + s-wave order parameter admixture, re- 
fiecting additional superconductivity on the chains. Since 
we expect the vortex nearest neighbors to lie along the 
nodal directions (e.g. [10]), the observed sign of distortion 
of the rhombic structure is qualitatively consistent with 
the phase-sensitive results [18]. With increasing field, we 
expect that the chain gap along b is further suppressed, 
so that the order parameter approaches (1^2 _y2, and the 
nodes move towards (110). Hence our observation of 



the vortex nearest neighbor directions approaching (110) 
with increasing field is evidence for a field-driven change 
in the nodal positions. This may be testable by high-field 
scanning tunneling microscopy measurements. 

In conclusion, we have used SANS to observe two field- 
driven first- order VL structure transitions in detwinned 
YBCO. We argue that the low-field transition at 2.0 T is 
driven by Fermi surface effects, and the high-field transi- 
tion at 6.7 T is driven by the order parameter anisotropy. 
The high-field transition is between hexagonal and rhom- 
bic structures with distortions of opposite sign, providing 
evidence for a d + s-wave admixture for the order param- 
eter. 
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ligen, Switzerland. We acknowledge financial support by 
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the SNE, and from the European Commission under 
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